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Abstract 

o : 
o 

^ An analysis based on the direct torque equations including tidal dissipation and a viscous 
^ core-mantle coupling is used to determine the damping time scales of O(IO^) years for free 
^ {^recession of the spin about the Cassini state and free libration in longitude for Mercury. 

The core-mantle coupling dominates the damping over the tides by one to two orders of 
^ tnagnitude for the plausible parameters chosen. The short damping times compared with 

the age of the solar system means we must find recent or on-going excitation mechanisms if 
^ Such free motions are found by the current radar experiments or the future measurement by 
C the MESSENGER and BepiColombo spacecraft that will orbit Mercury. We also show that 
^ the average precession rate is increased by about 30% over that obtained from the traditional 
O precession constant because of a spin-orbit resonance induced contribution by the C22 term 

in the expansion of the gravitational field. The C22 contribution also causes the path of 
^ the spin during the precession to be slightly elliptical with a variation in the precession rate 

that is a maximum when the obliquity is a minimum. An observable free precession will 
^ compromise the determination of obliquity of the Cassini state and hence of C/M^^iR^ for 

Mercury, but a detected free libration will not compromise the determination of the forced 

libration amplitude and thus the verification of a liquid core. 
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1. Introduction 

A major goal of the MESSENGER (Solomon, et al. 2001) and BepiColombo (Anselmi & 
Scoon, 2001) missions to Mercury and ground-based Radar Speckle Displacement Interferom- 
etry (RSDI) observations (Holin, 1988, 1992, 2003; Margot et al. 2003) is the determination 
of the structure and state of Mercury's interior. The probe of interior properties is based on 
the accurate determination of the four parameters, i, (f)Q, J2 and C22, where i is Mercury's 
obliquity, 0o is the amplitude of the physical libration in longitude, and J2 and C22 are the 
second degree harmonic coefficients in the expansion of Mercury's gravitational field. These 
parameters are sufficient for determining the factors in the equation (Peale, 1976a; 1981; 
1988; Peale et al. 2002) 
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where Mm and R are the mass and radius of Mercury, A < B < C are the principal moments 
of inertia, and Cm is the polar moment of inertia of the mantle alone. 

The first factor follows from the amplitude of the physical libration 

?, fB-A\ ^ . 959 4 \ 

where e is the orbital eccentricity, and Cm occurs in the denominator because a liquid core is 
not expected to follow the 88 day physical libration of the mantle. The second factor follows 
from 

and the third from (Peale, 1969) 
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MmR^ (sin I)/ic — cos / ' 

where ic is the obliquity of Cassini state 1, {E.g., Colombo, 1966:, Peale, 1969), which 
Mercury is expected to occupy {i — ic), I is the inclination of the orbit plane to the Laplace 
plane on which Mercury's orbit precesses at the constant rate — /x, and n is the orbital mean 
motion, and the equation assumes ic -C 1. In a Cassini state the spin axis is coplanar with 
the normal to the Laplace plane and the orbit normal as both the spin and the orbit normal 
precess around the Laplace plane normal with an approximately 280,000 year period. Fig. 1 
shows the configuration for Cassini state 1, where the ascending node of the equator plane on 
the orbit plane remains coincident with the ascending node of the orbit plane on the Laplace 
plane. In Fig. 1, the Z axis is normal to the orbit plane, the Y axis is in the orbit plane, and 
the X axis comes out of the paper. Eq. (4) shows the dependence of C/M^R'^ on J2, C22 
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and ic, and together, the last three equations show the necessity of knowing ic, 0o, J2 and 
C22 to determine C^/C. 

[Figure 1] 

Conditions for the success of the experiment are 

1) The hquid core must not follow the mantle during the 88 day forced libration in longitude. 

2) The core must follow the mantle during the 280,000 year orbit precession, where the spin 
axis is locked to this precession if it occupies the Cassini state. 

3) B — A must be due to the mantle alone. 

4) Adiabatic invariance must keep the spin close to the Cassini state during the slow variation 
of the orbital elements. 

The first two conditions are satisfied for viscous core-mantle couphng, topographical 
coupling due to core-mantle boundary (CMB) irregularities, magnetic core-mantle coupling, 
and gravitational coupling between an axially asymmetric solid inner core and the mantle 
(Peale, et al. 2002). The third condition requires that there not be long wavelength irregu- 
larities in the CMB such that dense core material contributes to B — A. Such irregularities 
could be induced by mantle convection (D. Stevenson, private communication, 2003), but a 
major contribution to B — A would occur only if the convective cells have circumferential 
scales comparable to the planet radius R. Otherwise, the plus and minus contributions due 
to short wavelength topography would tend to average to a zero total contribution to B — A. 
As the iron core must be near 0.75R in radius (Siegfried and Solomon, 1974), the horizontal 
scale of any convective cells is unlikely to greatly exceed the 0.25R thickness of the mantle. 
It has also been argued that Mercury's mantle convection ceased early in Mercury's history 
(Reese and Peterson, 2002), and the CMB should thereby be nearly axially symmetric with 
little or no contribution to B — A from the core. 

It is implicitly assumed in the above that all free motions associated with Mercury's spin 
are damped to near zero amplitude, and that the spin remains close to the current location 

of Cassini state 1 in spite of slow variation of the orbital parameters that define the state. 
The free motions include a free libration in longitude, a free precession of the spin about 
Cassini state 1, and a free wobble. The free libration in longitude results when the axis of 
minimum moment of inertia docs not point toward the Sun when Mercury is at perihelion. 
The gravitational torque on the permanent asymmetry of the planet, averaged around the 
orbit, tends to restore the alignment of the axis of minimum moment of inertia with the 
Sun at perihelion because of the 3:2 spin-orbit resonance, so viewed at perihelion, the long 
axis will tend to librate about the direction toward the Sun with a period of about 10 years. 
The free precession of the spin axis is like the 26,000 year precession of the Earth's spin axis 
about the ecliptic, except the precession is about the Cassini state instead of the orbit normal, 
where the deviation of the Cassini state from the orbit normal is less than 2 arcmin. The 
period of this precession is 0(1000 years). Free wobble or non principal axis rotation occurs 
when the spin axis does not coincide with the axis of maximum moment of inertia. The 
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spin axis then precesses around the axis of maximum moment of inertia in the body frame 
of reference with a period 0(500 years). The wobble is analogous to the Chandler wobble 
of the Earth. The angular momentum is conserved during free wobble, and the spin axis 
makes only small excursions relative to inertial space. The spin angular momentum is not 
conserved during the free libration and free precession as the planet responds to the external 
gravitational torque due to the Sun. In the former the spin axis direction is maintained as 
the magnitude of the spin varies periodically, whereas in the latter the spin describes a cone 
about the Cassini state, which is almost a cone relative to inertial space. 

The term "free" characterizes these motions as having arbitrary phases and amplitudes, 
which amplitudes, if significantly different from zero, could complicate the interpretation 
of the spacecraft and radar data constraining the properties of the interior. For example, 
a finite amplitude free precession would displace Mercury's spin vector from the Cassini 
state, and its arbitrary phase and amplitude would thereby increase the uncertainty in the 
obliquity of the state used in determining C/M^R^ in Eq. (4). Tidal friction has been shown 
to carry the spin to Cassini state 1 from almost any initial configuration on a time scale that 
is short compared to the age of the solar system (Peale, 1974; Ward, 1975). Free libration in 
longitude (9.2 year period, for {B — A) /Cm = 3.5 x 10~"^) is similarly damped, but such a free 
libration does not confuse the experiment, as the 88 day forced libration is superposed on the 
free libration (Peale, et al. 2002), and its amplitude is readily discernible. Any free wobble 
will also damp due to dissipation by internal friction as the equator plane changes relative 
to the body principal axis system in addition to tidal friction and core-mantle dissipation. 
For the Moon, tidal dissipation is about 4 times more effective than the dissipation caused 
by the changing orientation of the equator (Peale, 1976b). Finally it was asserted (Peale, 
1974) that the spin would remain close to the instantaneous position of the Cassini state if 
it were ever close because of the adiabatic invariance of an action integral that is equivalent 
to the solid angle swept out by the spin vector as it precesses around the Cassini state. 

If, on the other hand, free motions of measurable amplitude are found, excitation mech- 
anisms must be sought given the damping. Relatively small amplitudes of free motions can 
be detected by orbiting spacecraft (Zuber and Smith, 2000) and radar (Holin, 2003; Margot 
et al. 2003). The damping time scales constrain the types of excitation mechanisms that 
are likely to be responsible. For example, short time scales make excitation by the impact of 
large objects less likely (See Section 4). The existence of a liquid core introduces damping 
by core-mantle interaction that shortens the time scales for the damping of free motions 
deduced for tidal dissipation. This paper is the first of three where we consider the free 
motions and quantify the adiabatic invariance of the spin axis separation from the Cassini 
state. Here we derive and compare damping time scales for the two processes of tidal fric- 
tion and core-mantle dissipation with viscous coupling for free libration and free precession. 
Quantifying the adiabatic invariance and determining the damping of a free wobble require 
different approaches, and these will be considered in the following papers. The damping of 
the free rotational motions of the Moon have been determined earlier using a Hamiltonian 
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formulation with dissipation accounted for by introducing phase lags in the periodic terms 
describing the response of the Moon to various forcing potentials (Peale, 1976b). This analy- 
sis is general albeit somewhat complicated. For the libration and precession considered here 
we use a much simpler approach by fixing the orbit and writing the equations of motion 
directly in terms of the external and internal torques on the mantle and core. However, this 
simplified approach will not work well for the wobble, so the general Hamiltonian analysis 
will be used in a following paper to treat the wobble damping, although the torque equations 
will be applied there for the core-mantle coupling. In Sections 2 and 3, we show that the 
damping times for the free precession and the free libration in longitude are both near 10^ 
years with the core-mantle interaction being between one and two orders of magnitude more 
effective than tidal friction. In Section 4 we show that the time between impacts capable of 
generating observable amplitudes of free libration or precession exceeds 10^ years, making 
such excitation exceeding unlikely. Summary and conclusions follow in Section 5. We begin 
with the free precession, as the equations describing the damping of the libration in longitude 
are scalar versions of those describing the damping of the precession. 



2. Damping of free precession 

Generally, the spin axis precesses about Cassini state 1 which is offset approximately 
1.7 arcmin from the orbit normal. However, for the purposes of calculating the damping, we 
assume the orbit is fixed, so that the precession will be about the orbit normal. (The Cassini 
state for a fixed orbit can be thought to coincide with the orbit normal, since the state 
approaches the orbit normal as the orbital precession rate is decreased, or it can be thought 
not to exist at all.) Accordingly, we adopt an XYZ inertial system centered on Mercury 
with the XY plane coincident with the orbit plane, the X axis pointing from Mercury to 
the perihelion of the Sun's relative orbit, and with unit vector eo = e^ being the orbit 
normal. During the precession about the orbit normal, the damping of the precession angle 
about the orbit normal will be essentially the same as the damping of the precession about 
the Cassini state without involving the complications induced by the orbit precession. (See 
Peale, (1976b) for a complete analysis applied to the Moon.) 

We assume that the core mantle interaction is simply proportional to the difference in 
the vector spin angular velocities of the mantle and core as the simplest possible coupling. 
This coupling is consistent with a viscous coupling between two rigid spheres. Mercury is 
rotating slowly enough such that pressure coupling of the core and mantle from the elliptical 
shape of the core- mantle boundary is not important, which leaves viscous coupling the most 
probable. We also assume principal axis rotation. The equations describing the motions of 
the core and mantle are then 
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Cc^c = K^m-^c)^ (5) 

— * 

where ,^ are the angular velocities of the mantle and core respectively, C^.c are the mo- 
ments of inertia about the spin axis of the mantle and core, k is the coupling constant for 
the core-mantle torque, is the torque on Mercury's permanent figure due to the Sun, and 
T is the tidal torque. We first describe the free precession of Mercury, which is affected by 
the spin-orbit resonance {ipm) = 1.5n, where {ipm) is the average value of the magnitude of 
iprn, and n is the mean motion. 



2.1. Free Precession 

We start with the first of Eqs. (5), with T — k — and keep only the second order 
terms in the expansion of Mercury's gravitational potential Then the torque on Mercury 
is the negative of the torque on the Sun, so 



r X V$ = N = 



6C22 sin 9 sin 20 e^i 



-I- (3 J2 sin 9 cos 9 -\- 6C22 sin 9 cos 9 cos 20) e^ , (6) 

where G is the gravitational constant, and Mq are the masses of Mercury and the Sun 
respectively, R is Mercury's radius, J2 and C22 are second order harmonic coefficients in the 
expansion of Mercury's gravitational field, are the spherical polar coordinates of the Sun 
relative to the principal xyz axis system with ee and e<^ being unit vectors in the 9 and 
directions. The z axis is the axis of maximum moment of inertia coincident with the spin 
axis, and x and y are axes of minimum and intermediate moments of inertia in the equator 
plane. We define m = e^ as a unit vector along the spin axis of the mantle. The xyz axes 
are oriented with respect to the XYZ axes through Euler angles VL^i^ipm as shown in Fig. 
2, where Vt is the longitude of the ascending node of the equator on the orbit plane, i is the 
obliquity, and ipm is the angle between the node and the x axis in the equator plane, with 
the subscript m again referring to the mantle. 

[Figure 2] 

The expression for N in spherical coordinates relative to the body axes (Eq. (6)) is 
written in terms of the xyz components, and these are transformed to the XYZ components 
through the usual series of three rotations through the Euler angles. With the unit vector 
toward the Sun, f = cos /e^ + sin /e^ = sin 9 cos + sin 9 sin -|- cos 9 e^, where ex,Y,z 
and e.x^y^z are unit vectors along the respective axes, scalar products f ■ Gxi^yi^z yield the 
expressions for the products of circular functions of the spherical coordinates in N in terms 
of the Euler angles and the true anomaly / when ex^y^z are expressed in the XYZ system. 

— * 

These results and the transformation yield N in the XYZ system in terms of the Euler 
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angles and the true anomaly. 



N 



|3J2 smzcoszsm (/ — iZ)[sm/ — cos / e^J 

— 1.5J2sinisin2(/ — Q) 

+6C22sini[cosicos2'0sin (/ — Q) — cos(/ — Q) sin2'0](— sin/e^ + cos/e^) 

r 3 + cos 2i "I "1 

+6C22 ^ sin2(/ — r2) cos2-?/' + cosicos2(/ — f2) sin2-?/' e^j 



(7) 



To eliminate high frequency terms due to the orbital motions, we average Eq. (7) over an 
orbit period keeping everything constant except / and ijjm- With i/jm = l-5n {n — mean 
motion), we can write ip^ = V^mo + 1-5M, where M is the mean anomaly and ipmo is constant 
over the orbit period. With this substitution for ipm, we expand the circular functions in Eq. 
(7) to display products of circular functions of / and M with a^/r^. The terms yielding non 
zero averages over the orbit period are then 

1 dM 
53e3 



■df 



1 



2^3/2 



cos 3M ) = 



16 



+ 



^cos2/cos3M) = 



■ sin 2/ sin 3M 



123 

123 



-e e*" + 

2 16 



(8) 



where r = a(l — e^)/(l + ecos/) and dM/df = r^/ (a^Vl — e^) have been used. Some of 
these averages are easily done directly, whereas others are easily evaluated with the help of 
Cayley's (1859) tabulated expansions. 

When M = 0, the x axis is nearly aligned with the X axis as a condition of the spin-orbit 
resonance. This will be true for small i if -0^0 = With this substitution, the averaged 
torque on the permanent asymmetric figure of Mercury is 



(N) = 



3J2sinicosi 3C22 sini(l + cosi) 



2(1 -e2)3/2 
(cos Qex + sin Q e^) 



+ 



7 

-e 
2 



123 

f 

16 



X 



+ 



3C22sini(l + cosi) 53 3 , • o 
e (— cosize^ + smize^ 



16 

3C22 sini(l — cosi) 53 



16 



e (cos 3Q Bx + sin 3fl 



(9) 



where the Z components average to zero — the Z component with coefficient J2 because 
((a^/r^) cos 2/) = ((a^/r^) sin 2/) = 0, and the Z component with coefficient C22 because 
sin (2-?/'mo + 2fi) vanishes with ipmo = The first bracketed factor shows how the tradi- 
tional precession constant, involving only J2, has been modified by the spin-orbit resonance 
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with the addition of the C22 term, which decreases the precession period by nearly 30% from 
that due to the J2 term alone. (For the Earth, only J2 is important in determining the rate 
of precession, since terms involving C22 for the noncommensurate, rapid spin average to zero, 
and in addition, C22 ^ J2-) 

— * 

Before writing {N) in a form suitable for solving the equations for the spin vector motion, 
it is instructive to write Eq. (9) in a different way. The following identities are needed. 



cos VtGx — sin VLGy 

sin i (cos VLGx + sin fie^) 
sini(sinOex — cosOey) 



cos 2f2(cos VLBx + sin f2ey) 
+ sin 2fi(sin Vtex — cos Vt 
— m X 

-(m X ez) X Gz 



(10) 



If we rearrange terms and employ some trigonometric identities along with Eqs. (10) in Eq. 
(9), we can write 



— [Kicosi + Ks(l — cosi)]{ni X ez) 

—7^2 cos i cos 21] (m x e^) 
— i^2 sin 2f2(m x e^) x e^. 



The constants i^i,2,3 are defined by 
K2 = 



Sn^M^R^ „ 53 



2(1 - e2)3/2 + 3^^^ 1^2' 



123 
16' 



7 



C22 -e 

2C„Vrr V2 



123 , 



16 



(12) 



The m X terms in Eq. (11) describe the retrograde precession of m about the orbit 
normal, but now with a rate that varies with twice the precession frequency. The triple 
cross product term varies the obliquity i with the same frequency. Recall that the longitude 
of the ascending node of the equator on the orbit plane Q is measured from the perihelion 
direction. The precession rate is minimal when Q = Oorvr and maximal when Q = ±7r/2. 
The former corresponds to the node of the equator plane being parallel to the hue of apsides 
such that at the perihelion, the torque due the axial asymmetry of Mercury {K2 has C22 
as a factor.) vanishes. When Q = ±7r/2, the node is perpendicular to the line of apsides, 
and at perihelion, the long axis is tilted the maximum amount relative to the orbit plane. 
The torque and hence the averaged precession rate is thus maximal for this phase, since the 
biggest torque due to the egg shape of Mercury comes at the perihelion when the long axis 
is tipped above or below the plane of the orbit while still being oriented toward the Sun 
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as closely as possible because of the spin-orbit resonance. For Q — or n, the opportunity 
for the Sun to exert a large torque on the axial asymmetry at perihelion is lost, and so the 
average of the torque around the orbit is minimal. 

The triple product term in Eq. (11) lies along the projection of the spin vector ih on 
the orbit plane. So for 7r/2 > Q > the obliquity is increasing, whereas it is decreasing for 
> Q > — 7r/2. Hence, the obliquity has a minimum when the spin is crossing the X axis 
(Q = ±7r/2) and a maximum when it crosses the Y axis {Q = OorTr). The precession rate is 
maximal when the obliquity is minimal and vice versa, and the precession path of the spin 
is elliptical with the long axis along the Y axis. The reader will notice that there is also a 
component of the triple vector product torque that is parallel to m. This component may 
not be real because of the approximations involved in arriving at Eq. (9), but if it is real, 
there may be a displacement of the center of physical libration as a function of precessional 
phase. 

A good approximation to the precessional motion can be obtained by neglecting the 
last term in Eq. (9), since it is O(i^) smaller than the next smallest term. With rfi = 
sinisinQe^ — sinicosfiey + cosie^, we can write Eq. (9) in still another way with a 
rearrangement of terms. 



Kiruz + Ks{l - ruz) - K2 



[1 + nii 



rux ey. 



(13) 



where mx,Y,z are the respective components of ih. Then 



dm 



(N) 



where = 1.5n is fixed by the spin orbit resonance and can be removed from the time 
derivative. Prom Eq. (13) 



dlTLx 

dt 
dmy 
dt 



Kim, + KM - niz) - K', 



(l + "'z) 



Kini, + K^{l-niz) + K2 



2 

'1 + rUz 



ruy 



nix. 



(14) 



where for small obliquity, with ~ 1, the solution is 



niv — 



mPx cos 



Kl-KUt-rC 



K1 + K2 
Ki-K2 



K'^ 



(15) 
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where m° and ( are constants depending on initial conditions. With Cm — O-ISM^-R^, 
ip = 1.5n, J2 = 6x 10-^ and C22 = 1 x 10"^ (Anderson al. 1987), = 1.117 x lO-W/yr 
and K2 — 8.391 x 10~^rad/yr, the precession period of Mercury's spin about the orbit normal 

(or the Cassini state) is 562 years. If Cm = OAGM^R'^, the precession period would be 
reduced to 500 years. If the core were sohd and C = O.SSMmR'^ , the precession period would 
be 1031 years. This compares to a 1066 year period obtained numerically by Rambaux and 
Bois (2004). The ratio of the axes of the ellipse traced out by the projection of the spin 
vector on the orbit plane is m™"^/m™"^ = ^/(l{^~+l(2)J{Ki'-^^^ = 1.0075. This solution 
is consistent with the deductions about the precession from Eq. (11) discussed above. 



2.2. Tidal Torque 

The tidal torque is given by 



T= (f •f^)(f,xf), (16) 



j-6 

where k2 is the second degree potential Love number, f is a unit vector toward the Sun, and 
Tj. is a unit vector toward the tidal maximum, which is the sub-Sun position on Mercury a 
short time At in the past. 

= f - ^At, (17) 

where the time derivative is relative to the body system of coordinates. Replacement of fx 
with Eq. (17) yields 

r' 



3k2GM^R^At . , , 

r=^-f fxf, (18) 



where we have set f • = 1. So with the generic relation for a vector D. {dD / dt)body — 
{dD/dt) space — Co X D, and with D — > f = cos/e^ + sin/e^ and Co '0m(sinisinnex — 
sin i cos fley + cos i e^) + CIbz, it is easy to obtain 

fxf = ■0^sinicos(/ — f2)(— sin/ex + cos/ey) 

+(/ - Q - V'mCosi)e^ (19) 

Eq. (19) is substituted into Eq. (18) and averaged over the orbit period. Useful averages 
are 

P) = "(l + fe' + fe' + ^e»)/(l-r)« = n/:(e) 
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■' / 2 

/2(e) - /3(e) 



sinV = :ilw_^ow (20) 



where / = n(a^\/l — e'^)/r^ has been used in the first of Eqs. (20). The averaged tidal 
torque is thus 



+iPm Sm Z COS i I I ^ j 



e 



+ [^/i(e) - /2(e)^mCosi]e2|, (21) 

where we have neglected Cl compared to n and ipm- 

The relationship between the dissipation function Q and At follows from a simple ex- 
ample. The dissipation parameter Q for an system oscillating at frequency uj is defined by 
{e.g. Lambeck, 1980 p 14) 



dE , 
1 ——at 

1 r rft 



cuAt, (22) 



Q 2'kE* 

where the numerator is the energy dissipated during a complete cycle of oscillation and E* is 
the maximum energy stored during the cycle. For a tidally distorted, nearly spherical body 
with the disturbing body in a circular equatorial orbit, a cycle would consist of half a rotation 
of the distorted body relative to the body causing the tide. For a complex tide generated by 
a non circular, non equatorial orbit of the disturbing body, each periodic term in a Fourier 
decomposition of the tide would have its own maximum stored energy and dissipation over 
a complete cycle of oscillation. 

The response of an oscillator with forcing function F = A' sin Ljt when a; <^ cuoj with 
luq being the lowest frequency of free oscillation, can be represented by x — B' sin u!{t — At), 
where At is the phase lag in the response due to the dissipation as was assumed above for 
the tidal response. The rate at which the forcing function does work is dE/dt — Fx — 
A'B'u sin ujt cos uj{t - At) . Then 

^/ X f*^., A,^, [— COS (2LAjt — laJ At) sincuAt 1* 
E{t) = / Fxdt = A'B'u ^— ^ H 1 . (23) 



4a; 2 J ti 



The first term in Eq. (23) is the periodic storage of energy in the oscillator, and the second 
is the secular loss of energy. The maximum energy stored is just twice the coefficient of the 
cosine term, E* = A'B'/2, and the energy dissipated during a complete period of 27r/a; is 
AE — A'B'nujAt with sin a; At a; At. We use the energy stored as the energy increment 
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above the minimum energy in the first term since the stored tidal energy will always be an 
increase over the relaxed spherical shape of the body. Hence 1/Q — a; At as indicated in 
Eq. (22). Since is independent of frequency, Q is inversely proportional to frequency in 
this model. For Earth like materials, Q is almost independent of frequency for a wide range 
of frequencies above typical tidal frequencies (Knopoff, 1964), but to avoid a step function 
reversal in torque when a tidal frequency changes sign, it is not unlikely that the above 
frequency dependence will prevail for the small tidal frequencies. 

Q = 100 is typical of Earth-like materials (Knopoff, 1964) and apparently for Phobos 
generated tides on Mars (Smith and Born, 1976). For Mars the dominant frequency is 
uj = 2{il) — Hp), where ip is Mars rotational angular velocity and Up is Phobos' mean orbital 
motion. The factor 2 follows from there being two tidal cycles for each synodic period. With 
a rotation period of 24'*37'"23^ and an orbital period of 0.319 days, u = 3.14 x 10~^ rad/s 
leading to At 32 seconds for Q — 100. For Mercury, the fundamental tidal period is 
the orbit period, since Mercury rotates 180° relative to the Sun for each orbit. Although 
one might infer that Q is inversely proportional to frequency near a particular frequency to 
smooth the transition from positive to negative frequency, it cannot be the case that the 
same proportionality constant will apply over a very wide range of frequency as indicated 
by Q ~ 100 nearly independent of frequency between approximately 10 and 10^ cps for 
Earth materials (Knopoff, 1964). We can write At = P/{27rQo), with Qq being the value 
appropriate to the 88 day orbital period P. 



2.3. Damping the Obliquity 

To evaluate the damping of the free precession, or equivalently, the rate of change of 
the obliquity i of the mantle due to the torques on the right hand side of the first of Eqs. 
(5), we write cosi — Bg ■ ni. (Recall = e^.) With P = Cm4>m — C^mi^m^, 

— * — # — * 

If we note that d(3/dt = T, where T is the total torque acting on the mantle, and that 
dp/dt = m • r, Eq. (24) becomes 

— *_ — # 

dcosi -\-^ ^ . .e±-T 

-^-[eo- (e. • m)m] • - = smz-^, (25) 

where is a unit vector lying in the equator plane that is perpendicular to the node of the 
equator on the orbit plane. 



From Eqs. (11) and (25), 



sini — —K2sm icosisin2il, (26) 
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represents the periodic variation in the obhquity during the precession. Prom Eqs. (21) and 
with e± — — cos i sin Q + cos i cos Q + sin i e^, 



. - e^ • (T) 2 . 

Sim — = A4Sin I 



where 



nhie) - ^^^(/2(e) - /3(e) cos2Q) 



(27) 



= ? , (28) 

with 1/nQo replacing At as described above. The functions /j(e) are defined in Eqs. (20) 
with /i(e) = 1.72, /2(e) = 1.37, and /3(e) = 0.0779 for e = 0.206. 

We have only the core-mantle interaction term left to evaluate, which we determine with 
the analysis of Goldreich and Peale (1970). We note first that K2 <S Ki from Eqs. (12), and 

— * 

since the tidal torque and core-mantle interaction are small compared to A^, the zero-order 
solution of Eqs. (5) is approximately a uniform retrograde precession of the spin vector at 
the rate — O, where Vt is the magnitude of the motion of the projection of the spin vector on 
the orbit plane. This zero order motion of the mantle is determined from Eq. (13) 

—r- ^ i^i(rh X 60) = Kismi{—cosQex — sinQey), (29) 
where we have also omitted the Ks term because of the factor 1 — 0. Prom Eq. (29) 

271 

n = Ki = yr-\ (30) 

562 ' ^ ^ 

where the 562 year precession period is discussed above for = O-ISM^B?. If we substitute 
Eq. (30) into Eq. (29), restore the core- mantle interaction but continue to ignore the tidal 
torque, and divide both equations by C^fi^ to make them dimensionless, we can write 

— * 

where f = Qt, ■0m,c — '^m,cl^i ^' — ^/ (Cm^), and C — Cc/Cm- With S' — — ip^, we can 
write 

— * 

^ = i^'Jrh X e„) - k' S'. (32) 



If we assume t' — when D, — {i.e. with the ascending node of the equator plane 
on the orbit plane aligned with the perihelion direction), we can write ft — —Ctt — —t' and 



-14- 



replace Q, in the inertial system representation of ih x in Eq. (29) with —t'. Then 



d6'^ 
'dF 

dt' 



+ q;5^ = — ■0^sinicosf , 
+ aS'y — -^^sinisint', 



(33) 



where a — k'{C' + 1)/C" and 5'^ y,z the components of 5' along the respective coordinate 
axes. The solutions of these equations are 



5'^ — Ai exp —at' + tp'^ sin i 
S'y — A2 exp —at' + ip' sin i 



' —a cos t' — sin t' 

— cos f + n sill t' 
q;2 + 1 



(34) 



= ^3 exp — af , 

where are constants determined by initial conditions. 

With = — cosisinfie^ + cosicosi^ey + sinie^, and = — t', we substitute the 
steady state solution for 5' into the normalized form of Eq. (25), with 



k' sin i 



d'-e± 



k' sin^ i cos i 
q;2 + 1 '■ 



(35) 



to yield in dimensionlcss form 



d cos i f K2 /c' sin i cos i 

— ; — = sm I —— sm I cos i sm 2t H 

dt' \Ki a^ + 1 



+X4sin2i{n7i(e) - ^^^[/2(e) + /3(e) cos 2f]}, 



(36) 



where n' = n/Q. If we average Eq. (36) over a precession period, the circular functions with 
t' in their arguments vanish. Next with d cos i/dt' — — sin idi/dt' and setting sini w i and 
cosi ~ 1, we arrive at 



di 
dt' 



= —I 



(37) 



The time constant for the decay of the obliquity due to both a core-mantle interaction and 
tidal dissipation is the reciprocal of the coefficient of i in Eq. (37). 

We relate the coupling constant k' to the core viscosity by equating the time scale for 
the spin-up of a viscous liquid in a closed container with the time scale for relaxation of the 
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differential rotation. Prom Eqs. (34), tfie dimensioned time scale in the transient part of 
is 

1 C' Rc /„„s 

Ts=^ = ^ = —. — (38 

an {i + C')k'n (^„z/)V2' 

where the expression on the right is the time scale for a fluid with kinematic viscosity u, 
rotating differentially in a closed spherical container with radius Rc, to become synchronously 
rotating with the container at angular velocity ipm (Greenspan and Howard, 1963). Here 
tpjn — 1.5n, Rc — 0.75R is Mercury's core radius (Siegfried and Solomon, 1987), Q — 
(27r/562)yr-i and C" = CjCm ^ 1- So 

k' = = 8.59 X , (39) 

{i + C')nRc 

where u is expressed in cm^/sec. 

The kinematic viscosity of the Earth's core has been estimated from first principles to 
lie within the bounds 4.5 x 10~^ < u < 4.1 x 10"^cm^/sec by Wijs et al. (1998). If we 
choose V = 0.01 cm^/sec, then a = 2k' = 1.72 x 10^'^, and its square can be neglected in the 
denominator of the first term on the right of Eq. (37). With P = 88 days, R = 2439 km. 
Cm = O.ISMmR^, Mm = 3.302 x 10^^ g, the time scale for damping Mercury's free precession 
becomes 

89 years ^^^^ 



8.59 X 10- Vi/2 + 8.08 x 10-^k2/Qo ' 
Theoretical values of Mercury's k2 ~ 0.3 to 0.4 (Spohn et al. 2001; Van Hoolst and Jacobs, 
2003), and Qq ^ 100 is typical for Mars and Earth (Smith and Born, 1976; Knopoff, 1964) 
such that k2/Qo = 0.004 is a plausible value. With u = O.Olcm^/sec, we find Tp^ec ~ 
1.00 X 10^ years under the action of both tidal dissipation and core mantle interaction, and 
'Tprec — 1-04 X 10^ ycars and Tp^ec = 2.76 x 10^ years for core-mantle interaction and tidal 
dissipation respectively if they acted alone. 

The development of Eqs. (5) for numerical integration is given in Appendix A. An 
advantage of a numerical approach is that the consequences for the damping of large values 
of u, where the approximations used above are not valid, can be determined. Since, u ^ 0.01 
is unlikely, we shall not show the damping histories for large u, but remark that the damping 
time scale is a minimum for v lOOOcm^/sec where the precession amplitude is nearly 
critically damped. For larger u the couphng is sufficiently strong that the core nearly follows 
the mantle as both precess together and the precession period increases from 562 years 
to somewhat over 1000 years appropriate for a solid planet. The damping time scale for 
the precession amplitude is similarly increased from the minimum as u is increased. The 
Greenspan and Howard spin-up time scale in Eq. (38) also must be replaced by a viscous 
time scale for large u in relating u to k'. 



In Fig. 3 we demonstrate that the same time scale for damping the precession am- 
plitude from the combination of core-mantle and tidal dissipation is obtained numerically 
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as analytically for the same v — O.Olcm^/sec and k2/Qo — 0.004. If we set k — 0, the 
damping time estimated numerically similarly to Fig. 3 is 2.8 x 10^ years compared to 
2.76 X 10^ years determined analytically, which confirms the analytic development. Damping 
time scales for reasonable choices of u and k2 / Qo are all short compared to the solar system 
age, so we would expect the free precession to be completely damped, barring excitation by 
an unspecified mechanism. 

[Figure 3] 



3. Damping of free libration 

The treatment of the libration damping is considerably simpler than that of the preces- 
sion amplitude as we can assume zero obliquity as well as principal axis rotation. Eqs. (5) 
still apply, but are now scalar equations since there are only torques about the Z axis, which 
is now coincident with the z axis. With f x f = —{ipm ~ f)^z, we can write 

Cmi^m = N + T -k{7Pm-'>Pc), 

= k{ijm-ijc): (41) 

where C22 = {B — A) / {AMmR^) has been used. The relations between /, 0, and-j/'m are 
shown in Fig. 4, where the X axis is now drawn from the Sun to the perihelion with line 
SM drawn from the Sun to Mercury. The angles are hnked to the orbital motion through 
V'm = ^ + /; where / is the true anomaly. 

[Figure 4] 

Mercury's rotation deviates from iZ-v,,, = 1.5n by only a small amount, where i/;„t is 

the magnitude of -0^. We therefore write 7^ = ipm — l-5n as the deviation of the spin 
from the mean value. The angular deviation of, say. Mercury's axis of minimum moment 
of inertia ("long axis") from the position it would have had if the rotation were uniform is 
then 7^ = ijjm — 1.5M, where M = nt is the mean anomaly. 7^ is thus the angle between 
Mercury's long axis and the direction to the Sun when Mercury is at perihelion. The equation 
of rotational motion for the mantle becomes 

Ira + Sin (27„ + 3M - 2/) = - _ ^^), (42) 



where = GM^/a^ is used, with a being the orbit semimajor axis. The expansion of 
sin (27m + 3M — 2/) leads to factors (a^/r^) sin 2/ and [a? /r^) cos 2/, which when expanded 
in terms of the mean anomaly M with the help of Cayley's tables (1859) gives 

959 



B-A^.. 35- A 

7m + 3-- — G2oi(e)7,„ = -7^—7^ — 



l-lle^ + ^e^ + -- - \sint 
48 
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1 421 . , . ^ 
2e + ^e3 + ... Isin2t + 

+ 7^-7^(^-^), (43) 

where e is the orbital eccentricity, G2Qi{e) — 7e/2 — 123e^/16 H is a Kaula (1966) eccen- 
tricity function. We have made the equation dimensionless by dividing by and letting 
nt = M ^ t. Also, sin27TO — >■ 27^, cos27m, 1 have been used, and a Fourier series in t 
multiplied by 7^ has been omitted from the right hand side of the equation, since it is so 
much smaller than the unit coefficient of the displayed series. 

If we neglect the weak tidal and corc-mantlc torques for the moment, Eq.(43) is a simple 
forced harmonic oscillator equation describing the free and forced librations of longitude for 
Mercury. The amplitude of the sin2t term is about 11% of that of the sint term and 
higher order terms arc even less. Since the overall amplitude of the physical libration is not 
changed by this second term in the series, and higher order terms contribute negligibly, we 
retain only the sint in defining the amplitude of the physical libration in Eq. (2). Although 
this truncation of the series gives a reasonably good approximation to the physical libration 
angle, it gives a poor representation of the deviation of the angular velocity from the mean 
value. This is shown in Fig. 5 where the flat tops on the angular velocity curve at perihelion 
result from the fact that the orbital angular velocity is close to the spin angular velocity 
near perihelion, and with the long axis nearly aligned with the direction to the Sun during 
this time, the torque on the axial asymmetry is nearly zero. A detailed explanation of the 
behavior of dip^/dt near the perihelion is given in Appendix B. 

[Figure 5] 

The free hbration is treated by averaging Eq. (43) over an orbit period such that the 
periodic terms responsible for the physical libration vanish. T is replaced by its averaged 
value (T) and 7^ and 7c arc uniform on the orbital time scale. We have already specified the 
form of the frequency dependent tidal torque in Eqs. (41), where At (defined in Eq. (17)) 
is expressed in terms of ipm — f ■ 



1 r 



2 n / 



2t:GQ,VI^ Jo r 

-f(v + -^), (44) 



n ^ 

where we have used the second form of the mantle equation in Eqs. (41), dM/df — 
I (a^Vl — e^), / = no? yl — e^ / r"^ , and where 
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^ 3 (1 + 15(:V2+ 15^78 + 5^71(3) 3 /i(e) 



2 (l-e2)3/2(i + 3e2 + 3e4/8) 2 /2(e)' 

which agrees with Eq. (24) of Goldreich and Peale (1966), except for the correction of the 
coefficient of from 3/16 to 5/16. 

We divide Eq. (43) by = 3{B — A)G2oi{e) / Cm and further normalize time by xt = t'. 
Time t — > t' is now normalized by the free libration frequency instead of the precession 
frequency Q as in Section 2. The average of Eq. (43) then becomes 



" Im 
dt'^ 

dH k' 



k' 

^^,2 = ^(7m-7;), (45) 



where F' = F/CmX^n'^, k' = k/CmXn, C = Cc/Cm and Ymc — im,c/xn. Similar to the 
treatment of the precession, the fact that the tidal and core-mantle torques on the mantle 
are small means that the mantle will librate almost as it would if k and F were zero. We 
follow the same procedure as in Section 2 but now in scalar form. The zero order solution of 
the first of Eqs. (45) is then 7^ = TmoSinf and = 7;„ocosf, where (7m, 7m) = (0,7mo) 
when t' = 0. Subtracting the second of Eqs. (45) from the first yields an equation for 

<^ = 7m - 

§ + k'^^S = -y^, sin t' -F'{V + xi^, cos t') , (46) 

where we have substituted the zero order solutions for 7^ and 7^ on the right. Now F' — 
1.5 X 10-2A;2/Qo, where we have used x = ^3(5 - A){le/2 - 123eVl6)/C^ = 0.0262, 
Cm = 0.18MMi^^ inferred from Siegfried and Solomon (1974), {B - A) /Cm = 3.5 x 10"^, 
and Mercury's mass Mm = 3.302 x 10^^ g. The coefficient of cost' in Eq. (46) is a factor of 
3.9 X 10~^ less than the coefficient of sinf so we can neglect the former term. We keep the 
constant term leading to the solution 

8 = A4exp (—at) + 7?no ^ asint) ^, (47) 

1 + a 

where a — k'{l + C')/C' and A4 is an arbitrary constant. 

The steady state solution for 6 is substituted back into the first of Eqs. (45) to yield 

''^?+(T^ + f'-)%+ (1-^^)7;^-^. (48) 



dt'^ \l + a^ J dt' \ 1 + aV ^ + 

where we have used the zero order solutions, 7^ = 7^0 ^i^^' ^^<i 7m — i'mo^^^^' 
steady state solution for S and rearranged terms. Tidal friction displaces the zero point of 
the libration by -F'V/{1 + C') = 1.83 x lO-^fca/Qo = 7.34 x 10"^ radians or ~ 1.5 arcsec for 
^2/Qo = 0.004 {x = 0.0262, V = 0.244). A measure of this displacement would determine 
k2/Qo for Mercury. (Tidal friction attempts to slow Mercury's spin to a value that is less 
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than the current 1.5n, so at perihelion, the axis of minimum moment of inertia is displaced 
from exact alignment with the Sun (7^ = 0) until the averaged torque on the permanent 
deformation balances the tidal torque and keeps the spin at the resonant rate.) 

The time scale Tub for damping the free libration amplitude is twice the inverse of the 
coefficient of 7^ in Eq. (48). The transient part of the expression for 5 in Eq. (47) is 
the same as that for the precession analysis, so the coupling constant k' is related to the 
kinematic viscosity of the core by Eq. (39), except Cl must be replaced by xn. Then 

The dimensioned time scale for damping the free libration is given by 

2 

Tub 



{k'/{l + a'^) +F'x)xn 
2 



{k' + F'x)xn 

2.92 rmyears 
(1.41 X 10-W2 + 3.93 X 10-^k2/Qo 



(50) 



from Eq. (48), where a = 2k' is small and its square can be neglected. Wc again choose 
u = 0.01 cm^/sec after Wijs et al. (1998) along with k^jQ^ = 0.004 discussed above, we find 
Tub ~ 1-86 X 10^ years under the action of both tidal dissipation and core mantle interaction, 
and Tub — 2.07 x 10^ years and Tub = 1-86 x 10^ years for core-mantle interaction and tidal 
dissipation respectively if they acted alone. Eqs. (41) are numerically integrated in Fig. 6 
for u — O.Olcm^/sec and k2/Qo — 0.004, where the short period fluctuations in Fig. 6 are 
the 88 day forced librations superposed on the free libration. The approximately determined 
time scale is 1.8 x 10^ years, which is close to the analytic value of 1.86 x 10^ years. Damping 
time scales for reasonable choices of u and k2/Qo are all short compared to the solar system 
age, so we would expect the free librations to be completely damped, barring excitation by 
an unspecified mechanism. 

[Figure 6] 



4. Consequences of Impact 

This is an exercise to determine what size impactor would be necessary to generate an 
observable free libration or free precession. We shall assume Mercury has no free motions 
prior to the impact, and that there is a complete transfer to Mercury of the entire angular 
momentum of the impactor relative to the center of mass of Mercury. If the impact ejecta 
is distributed symmetrically about the point of impact, it will carry away little angular 
momentum. This approximate conservation of angular momentum allows us to infer that 
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the angular increment in the amplitude of the precession is simply the component of the 
impactor angular momentum that is perpendicular to the spin axis divided by the spin 
angular momentum C^iprn- The increment in the free libration amplitude follows easily 
from the parallel component of the impactor angular momentum. Wc shall assume that a 
free libration amplitude or free precession amplitude of 0.1 arcmin is detectable. 

Generally, an impactor, which can be either a comet or asteroid, will strike at a point 
with spherical coordinates R, 6, (j) relative to the body principal axis system with velocity 
relative to the center of mass of Mercury oiV = 1/^6^ + 1^6^ + V^e^. The angular momentum 
of the impactor relative to the center of mass of Mercury is given by 

L,^RxMV = MR{Vge^ - V^ee), (51) 

where Ai is the mass of the impactor. With e<^ = — sin 4> e-^+cos (f) By and = cos 9 cos (f) 6-^+ 
cos 9 sin (pBy — sin ^ e^, we can write 

Li — M.R[{— Vq sin (f) — V^i, cos 9 cos (l))ex 
+{Vg COS — COS 9 sin (f))ey 

+V^sin9e,], (52) 

— * — * — * 

With the increment in Mercury's angular momentum, AL^^ = Lj, the components of L^ = 
L^j^j + ALm along and e^^ indicates that a free wobble has been induced. But that is not all, 
since the redistribution of mass due to the crater and its ejecta redefines the inertia tensor, 
so x, z arc no longer the principal axes. The new axis of maximum moment being offset 
from the original z axis also contributes to the wobble. The component AL^ means ip^n has 
been incremented thereby creating a free libration in longitude. Transforming L^ to the 
inertial system defines a new direction of the total spin angular momentum which no longer 
coincides with the Cassini state, which has itself been changed because of the change in the 

— * 

inertia tensor. Lm now precesses around the new position of the Cassini state. A random 
impact thus excites all three free motions (see Peale, 1975 for a more detailed treatment), 
but how large must the impact be for detectable amplitudes? 

First we note that 

Lj|| = MRV^ sin 9 

L,^ = MRiV^ + cos" 9f/^ (53) 

are the components of impactor angular momentum just before impact at coordinates {9, 0) 
that are parallel and perpendicular to the spin axis. We consider first the excitation of 
a free libration, which involves only L^y. The free libration is characterized by Eq. (43), 
where we neglect the right hand side to yield a small amplitude pendulum equation, whose 
frequency of oscillation is xn = n^J'i^B — A)G2m{e) /Cm- Recall that 7^ = V'm ~ 1.5M is 
the angular deviation of the x axis from the position it would have had if Mercury were 
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rotating uniformly at angular velocity 1.5n. If 7^ = and 7^ = 7^ at i = 0, the amplitude 
of the libration is 

„,Tnax -,0 /, , "^^11 l'KA\ 

Im = = -^T— > (54) 

where the final form assumes that all of the libration amplitude is induced by the impact. 

If we set 7^"^ = 0.1' = 2.91 x 10~^rad as the minimum observable, choose V^p = 
60 km/sec as a common impact velocity, set Cm = 0.18Mm-R^ as before, choose (B—A) / Cm = 
3.5 X 10^^ so that xn = 0.026n, and maximize the effect of the collision for a given Ai by 
assuming 9 = tt/2, we find the necessary A4 = 1.51 x 10^^ g. For a density of 1 g/cm^, the 
radius of the impactor is 0.712 km. 

Holsapple (1993) has described a careful development of scaling laws that predict the 
properties of impact craters from the size and velocity of the impactor and the strength 
and gravitational acceleration of the body on which the impact occurs. His Eq. (22b) gives 
the rim radius for large craters where gravity dominates the strength of the materials in 
controlling the outcome of the crater formation. 

Rr = 10.14G"-°-i^a'°-«3f;0.34^ ^55^ 

where Rr is the rim radius of the crater just after it is formed in meters, G' is the ratio 
of the gravitational acceleration at the surface of the body to that of the Earth, a' is the 
radius of the impactor in meters, and U is the velocity in km/sec. Eq. (55) ignores the 
ratio of the impactor density and surface material density raised to a small fractional power. 
With Mercury's gravitational acceleration of 363 cm/sec compared to Earth's 980 cm/sec, 
a' = 712 m and U = = 60 km/sec, we find Rr=11.3 km. For a crater this large, it gets 
even larger because of slumping of the walls, with the expression for the final crater size in 
terms of Rr given by Eq. (28a) of Holsapple (1993). This is called a transition from a simple 
bowl-shaped crater to a complex crater with a shallow floor and central peak. 

f = if) ' 

where Rf is the final rim radius, and i?* is the transition radius from simple to complex 
craters. R^ is a function of the surface gravity with i?* = 5 km for Mercury (Pike, 1988) (9.5 
km for the Moon and 1.5 km for Earth and Venus). Substitution of this value of R^ into 
Eq. (56) yields 12.4 km as the rim radius expected for the impact event leading to a 0.1' 
free libration amplitude for Mercury. This is a minimum value because of the choice of the 
extreme optimum values of and sin 6* in Eq. (54). It should be noted that the Holsapple 
crater radii are for vertical impacts, whereas we have calculated the angular momentum 
impulse for an almost grazing impact. However, this caveat does not change the conclusion 
that even the least detectable amplitude of a free hbration of 0.1' excited by an impact would 
leave a more than 20 km diameter pristine scar on the surface of Mercury. 
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For an impact excitation of the precession, the angular amphtude would be 

Ai = ^ = (57) 

where we have assumed the impact is along a meridian for simplicity. If Ai = 0.1', Vg = 
60 km/sec, and Cm = 0.18Mm-R^, then M = 8.72 x 10^^ g. The impactor radius is then 
a' = 2.75 km for a density of Ig/cm^. It takes a much more energetic impact to excite an 
observable free precession than an observable free libration. These numbers, when substi- 
tuted into Eqs. (55) and (56), yield Rf = 41.7 km, and the pristine crater left behind from an 
impact causing a 0.1' increment in the amplitude of the free precession would be more than 
80 km in diameter. The fact that the Cassini state itself has shifted shghtly because of the 
redefinition of the inertia tensor would not change the free precession amplitude sufficiently 
to alter this conclusion. 

In Table 1 of Levison and Duncan (1997), the rate of impact of Jupiter family comets 
on Mercury is given as 8.6 x 10~^° comets/year. These are for comets of absolute magnitude 
< 9, which correspond to radii a' > 1km (H. Levison, private communication 2005). The 
collision rate of very long period comets is probably small compared to this rate for the 
Jupiter family comets. (H. Levison, private communication 2005). More than a billion 
years would elapse on average between impacts on Mercury of sufficient size to excite a 0.1' 
amplitude of the free libration in longitude, and because the larger impactor required, even 
more time between excitations of a 0.1' amplitude of the free precession. Comparing this 
time between excitations with the O(IO^) year damping times for the libration and precession 
amplitudes means collisions almost certainly could not be the cause of any observable free 
libration or free precession. 



5. Summciry and Conclusions 

The spin orbit resonance of Mercury leads to a 30% reduction of the period of precession 
of the spin about the Cassini state from a contribution of the C22 term, compared to the 
usually derived precession period depending on J2 alone. This effect also leads to an elliptical 
precessional path of the spin, where the ratio of the axes of the eUipse is about 1.0075. 
There is a slight variation in the precession rate about the elliptical path with the maximum 
occurring when the separation of the spin from the Cassini state is minimal. 

The simultaneous effects of tidal dissipation and a viscous liquid core-solid mantle cou- 
pling lead to time scales for damping Mercury's free libration in longitude and free precession 
of the spin about Cassini state 1 that are near 10^ years. Because Mercury's liquid core is so 
large, the time scale for damping the precession amplitude by the core-mantle coupling acting 
alone is a factor of about 25 smaller than the time scale due to the tides alone for plausible 
choices of the parameters. For the libration damping, the core-mantle time scale is a factor 
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of about 10 smaller than that for tidal friction. For all reasonable choices of the parameters, 
the time scales are short compared with the age of the solar system. So if a free libration 
or free precession amplitude is found by the radar experiments or by the observations of 
the MESSENGER or BepiColombo spacecraft, still unspecified recent or ongoing excitation 
mechanisms must be sought. It is possible that there exist observable forced deviations of 
the spin axis from the Cassini state. If the forcing terms are identified, the amplitudes and 
phases of such deviations should be known and should not hinder the identification of the 
Cassini state obliquity. On the other hand, the forcing terms are likely to be of fairly long 
period, as is hkely if the amplitudes of the deviations from the Cassini state are observable. 
As such their identity may not be obvious and the Cassini state position could be thereby 
uncertain until careful numerical integrations mapped such forced motions. 

An observable libration amplitude (assumed to be 0.1') is much easier to generate by 
impact than an observable precession. But such impact excitation of either motion is very 
unlikely, since there is an estimated average time span between impacts of sufficient size of 
about 10^ years with current cometary fluxes. Excitation by impact no more than a few 
damping time scales in the past would leave a fresh crater larger than 20 km diameter for 
the excitation of an observable amplitude of libration and larger than 80 km diameter for the 
excitation of an observable precession amplitude. An observable free precession will make 
the obliquity of the Cassini state, and thereby the determination of C/MmR^ uncertain. 
But an observable free libration does not hamper the determination of the physical libration 
amphtude. If Mercury's core were solid, the term with the kinematic viscosity would vanish 
in the expressions for the damping times of the precession and libration (Eqs. (40) and (50). 
But the resulting damping times due to the tides alone would be longer that these equations 
would indicate for i/ = 0, because the constants would then involve the moment of inertia of 
the entire planet instead of just that of the mantle. 

This analysis apphes to other bodies in spin-orbit resonances such as the Moon, the 
regular satellites of the major planets, and Pluto's satellite Charon. The terms selected 
from the expansion of the torque expressions by the averaging process would correspond 
to rotation synchronous with the the orbital mean motions instead of 1.5 times the mean 
motion used here. There are similar contributions to the rate of precession about the Cassini 
state from C22, since the synchronous rotation ensures terms involving C22 do not average 
to zero and C22 will be comparable to J2 in magnitude. There are indications that some 
of these bodies even have liquid layers in their interiors {e.g., Spohn and Schubert, 2003), 
although Europa would be a special case with a liquid ocean [e.g. Greenberg, 2004) 

It will be more fun to seek alternative dynamic mechanisms of excitation of observed free 
motions of Mercury than if no such motions are found. A possible albeit unlikely excitation 
of libration may lie in a coupling between the torques on the permanent deformation as 
affected by the spin-orbit resonance and the small obliquity of the Cassini state. Although 
such a libration would not be "free," a near resonant forcing frequency could mimic a free 
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libration. The investigation of such speculation is a future objective. 



6. Appendix A: Numerical Integration of the Precession Damping 

We want to integrate Eqs. (5) over the precession time scale without following high 
frequency motions on the orbit or libration time scales. The averaged torque on the perma- 
nent deformation (N) is given in Eq. (9), where we set V'm = + 1-5M with t/jQ — — Q to 
accommodate the conditions of the spin-orbit resonance in the averaging process. The last 
term of Eq. (9) is of order i'^ smaller than the next smallest term, and it is neglected in Eq. 
(13) written in terms of the components of the unit vector along the spin axis ih. It is this 
latter form of the averaged torque on the permanent deformation that we will use here. 

— * 

The averaged tidal torque is given by Eq. (21). The component of (T) parallel to the 
spin axis is balanced by a slight offset of the axis of minimum moment of inertia from the 

— * 

direction to the Sun at perihelion. However, we have removed that degree of freedom in (N) 
by assigning ■0o = after Eq. (9). In order to avoid a secular change in i/jm, we keep 
only the component of (T) perpendicular both to the spin vector and to the line of nodes. 
The unit vector in this direction is defined in Eq. (25). From this last equation, we see 
that this is the only component of (T) that affects the obliquity. With ipQ = —fl and the 
elimination of the tidal acceleration of the spin, we have eliminated integration over the free 
libration period leaving only the precession time scale in the calculation. So in place of (T), 
we shall use 



-ex = -fV4 sm I 



nh{e) - ^^^^{h{e)-h{e)cos2Q) 



X 



(— COS i sin Vtex + cos i cos ^2 e^ -|- sin i e^) (58) 

where is defined by Eq. (28) with l/{nQo) replacing At as described in Section 2.2. 
Inclusion of the /3(e) part of the bracketed coefficient in Eq. (58) leads to a sini in the 
denominator, which complicates the calculation. However, /3(e) ~ 0.056/2(e) for e = 0.206 
so its omission will lead to a small error in the motion of the spin axis and an indiscernible 
error in the damping time scale. With this additional approximation, we arrive at 

(T) • ex 

• "" e^ = K4[nfi{e) - V'm"^z/2(e)][-mzmx e^ - mzmYey + {ml + ml) e^], (59) 

where the components of the unit vector along the spin axis of the mantle, m. 

With d{Cra'4'm)/dt = Cm'ipmdm/dt where ipm = 3n/2 = constant, we can average Eq. (5) 
over an orbit period, insert the averaged expressions (N) and the component of (T) parallel 
ex from Eqs. (13) and (59), divide by Ki to make the equations dimensionless, set Kit — t' 
and k/{CmKi) — k', and write six scalar equations in six unknowns, mx,Y,z,Cx,Y,z, where 
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Cx,Y,z are the components of the unit vector parallel to the spin of the core. 
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IF 



drUz 
~dF 
dcx 
dt' 
dcY 

HF 
dcz 
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■^{mx-Cx), 
^(my-Cy), 
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mzmy — k'{mY — Cy), 
{m\ + ml) - k'{mz - Cz), 



(60) 



are the equations that are integrated numerically and that describe the damping of the 
precession amplitude due to a core-mantle interaction and tidal dissipation. Here C — 
Cc/Cjn, and we have set ipm — l-5n in the term from the tides. 



7. Appendix B: Description of torque reversal near perihelion. 

This is an explanation of the nearly flat top on d'y^/dt in Fig. 5. In this figure there is 
a local minimum in 7^ at perihelion and two local maxima at equal times before and after 
perihelion. The angular acceleration is zero at these extremes and corresponds to the axis of 
minimum moment of inertia (long axis) of Mercury being aligned with the direction to the 
Sun. We ignore the slight variation in rotation due to the physical librations in determining 
the torque on Mercury by assuming it is rotating uniformly with tjjm = 1.5n. At perihelion, 
the long axis is also pointing toward the Sun, but / > ipm = l-5n at this point in the orbit. 
The long axis does not keep up with the motion of the Sun in Mercury's sky as Mercury 
passes perihelion, so it starts to point away from the Sun in the direction of the orbital 
motion. The angle between the Sun and the long axis continues to grow and reaches a local 
maximum when / = which occurs when / = 25.743° and M = 16.777°, 4.101 days after 
perihelion passage. The torque on Mercury during this time tends to increase the angular 
velocity corresponding to the positive slope on d'ym/dt in Fig. 5 just after perihelion passage. 
The increasing slope has an inflection at and starts to decrease after the time when / = ijjm, 
which is when the angle between the long axis and the Sun starts to decrease. The slope 
of the dipm/dt — d'jm/dt curve reaches zero when the axes are again aligned. To bring the 



-26- 



long axis back into alignment with the Sun after perihehon passage, Mercury must rotate 
through the same angle as it has moved in its orbit, so that / = 1.5M at the time of the 
first post perihelion alignment. This corresponds to / = 44.442°, M = 29.628° at 7.242 days 
after perihelion passage. After this time, the long axis points on the other side of the Sun 
and the angular velocity starts to decrease. 

The motion is symmetric with respect to the perihelion passage so the numbers corre- 
sponding to the approach to perihelion are the negatives of those above. It is ascertained 
that Mercury is almost synchronously rotating with its instantaneous orbital motion near 
the perihelion, during which time the long axis never deviates much from pointing toward 
the Sun. Hence, 7^ is almost constant for a time near the perihelion as shown in Fig. 5, 
with the small deviations from a constant value explained above. 
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Fig. 1. — Geometry of Cassini state 1. The Z axis is perpendicular to the orbit plane, 
and the ascending node of the orbit on the Laplace plane and the ascending node of the 
equator on the orbit plane both coincide with the X axis, which comes out of the paper. 
The obliquity of the Cassini state is ic- 
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Fig. 2. — Coordinates for describing Mercury's precession. The XY plane is the orbit plane 
with the X axis pointing toward the perihelion position, Vl is the longitude of the ascending 
node of the equator plane on the orbit plane, i is the obliquity, and ipm locates the x body 
axis of minimum moment of inertia (of the mantle) relative to the ascending node. The z 
axis is the axis of maximum moment of inertia, which is also the spin axis. Spherical polar 
coordinates rOcj) locate the Sun relative to the body axes. 



CO 

. I— I 

a 
o 

W 
-i-J 

0) 

c 
o 
a 

B 
o 
o 



0.02 



-0.02 - 



1 1 1 1 1 1 1 1 

A / 


1 1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
Vec=^t/in(A/yi; - 






'\ 


- 


l.OxlO^ yrs- 










- A 






.4 /- 






u = 0.01 cm 


Vsec — 


1 1 1 1 1 1 1 1 


/ v : 

1 1 1 1 1 1 1 1 1 1 1 1 1 


k/Q= 0.004 



0.5 1 1.5 2 2.985 2.99 2.995 3 
t (Precession periods) txlO"^ 



Fig. 3. — Numerically determined damping of Mercury's free precession with core kinematic 
viscosity u = O.Olcm^/sec and tidal factor k2/Qo = 0.004. rrix and Cx label the curves for 
the X components of the unit vectors along the mantle and core spin axes respectively. 
is a small but otherwise arbitrary initial amplitude. 
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Fig. 4. — Angles used in the discussion of libration in longitude. SX is directed from the 
Sun to the perihelion of Mercury's orbit and Mm denotes Mercury, where x is the axis of 
minimum moment of inertia. 
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Fig. 5. — Forced libration of Mercury showing the deviations of the motions from pure 
harmonic form. The fiat portion of d^m/ dt results from the fact that the rotation is nearly 
synchronous with the orbital motion for a range of true anomaly / on either side of the 
perihelion while the axis of minimum moment of inertia is nearly pointing to the Sun. (See 
Appendix B.) 
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Fig. 6. — Numerically determined damping of Mercury's free libration in longitude with core 
kinematic viscosity v — O.Olcm^/sec and tidal factor k2/Qo = 0.004. The high frequency 
oscillations superposed on the free libration are the 88 day forced hbrations. {B — A)/ Cjn — 
3.5 X 10~^ determines the free hbration period of 9.2 years. Aq is a small but otherwise 
arbitrary initial amplitude. 



